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The high energy spectral shape of 3C 273 is usually understood in terms of Inverse-Compton 
emission in a relativistic leptonic jet. This model predicts variability patterns and delays which 
could be tested if simultaneous observations are available from the infrared to the GeV range. 
The instruments IBIS, SPI, JEM-X on board INTEGRAL, PCA on boai'd RXTE and LAT on 
board Fermi have enough sensitivity to follow the spectral variability from the keV to the GeV 
and to compare them with model predictions. We are presenting preliminary results on the high 
energy spectrum of 3C 273 and its variabihty and compare these results to predictions. We found 
that a single component is not able to adequately fit the multiwavelength spectrum (5 keV - 10 
GeV) when the statistics in Fermi data is sufficient to constraint 7-ray emission. This suggests 
that the X-rays emission is not originated in the jet. A possible explanation could be that X-rays 
are dominated by a Seyfert-like thermal inverse Compton. 
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Figure 1: Top panel: PCA lightcurve in 0.2 - 75 keV. Bottom panel: Fermi lightcurve in 0.1 - 100 GeV. 
Area marked by dashed lines shows the selected flaring time epochs in both lightcurves. 

1. Introduction 

The bright quasar 3C 273 is one of the best monitored Active Galactic Nuclei (AGN). It shows 
several interesting blazar-like features (like the presence of a jet with superluminal motion and high 
variability) and Seyfert-like features (like a blue bump and variable emission lines). Many works 
have been published to describe in details its spectral and temporal characteristics, and a wide 
database of multiwavelength observations has been collected and it is available online. ' The aim of 
this work is to look for evidence of spectral variability related to different states of the brightness 
of the source. We analysed the high energy spectrum, in the keV to GeV energy range, of 3C 273 
using data from RXTE-PCA (2 - 70 keV), JEM-X, IBIS/ISGRI and SPI on board INTEGRAL (3 
keV - 1 MeV) and Fermi-LAT (100 MeV - 10 GeV). We selected two different epochs based on 
the behaviour of PCA and LAT lightcurves: the first epoch corresponds to the sum of the several 
flaring times in PCA while the second to the big flares in the LAT lightcurve (Fig. 1). INTEGRAL 
data were accumulated on the PCA and LAT flaring epochs to build the multiwavelength spectra 
of the source. 



2. Data Sample and Selection 

PCA spectra and lightcurves have been taken from the HEAVENS web interface.^ The signal- 
to-noise ratio of PCA is good enough to detect a well structured lightcurves, with clear flaring 
structures (Fig. 1 top panel). We selected the time intervals in which the source is flaring, defining 
the PCA flaring epoch, i.e. when the source is flaring in X-rays. 

Similarly we selected the LAT flaring epoch, i.e. when the source is flaring in /-rays, looking 
at the Fermi lightcurve (Fig. 1 bottom panel). The LAT lightcurve is indeed different from the 
PCA lightcurve: it shows a single big flare with some structures followed by another shorter flare. 
The first flare is described in detail in [1]. Except for this outburst the count rate is relatively stable 
showing few structures and few isolated points with high count rate. 



'http://www.isdc.unige.ch/3c273/ 
^http://www.isdc.unige.ch/heavens/ 
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The INTEGRAL lightcurves do not exhibit structures Uke the PCA Ughtcurves, so they are not 
used to select intervals of interest. INTEGRAL data have been used for spectroscopical purpose 
only. 

For both the two epochs described above a multiwavelength spectrum with data from all the 
instruments (PCA, JEM-X, IBIS/ISGRI, SPI and LAT) has been buih. 



3. Spectral Analysis 

We fitted the data with several models and compared the results. We used Broken power- 
law, Cutoff powerlaw, Log-parabola and Cutoff powerlaw-i-Log-parabola model, where the Log- 
parabola spectral law is 



N{E)=N,/ 



IkeV 

being Nip the photon flux at energy E = I keV. 

The multiwavelength spectral energy distribution (SED) with the best fit models for Broken 
powerlaw and Cutoff powerlaw+Log-parabola are shown in Fig. 2. In both panels, blue points 
and line denote the spectrum at LAT flaring time (henceforth spectrum A), yellow points and line 
at PCA flaring time (henceforth spectrum B) while the green points data taken from Litchi et al. 
1995 [2] are related to a multiwavelength campaign on 3C 273 performed in 1991, and they are 
used here for comparison (henceforth spectrum C). The observations in the X and y-ray bands were 
performed by GINGA and OSSE, COMPTEL and EGRET on board CGRO. Coloured solid lines 
are the respective best fit models and dashed lines are the individual additive components of the 
model (see section 3.2). A summary of the best fit parameters for the models is reported in Table 
1. 

Cutoff powerlaw model is not discussed in detail because it cleary does not fit the data. Fermi 
data can not be modeled by an exponential cutoff. For the Cutoff powerlaw model we obtain 
^red — 4^ f^'^ spectrum A and xjed — 37 for spectrum B. 

3.1 Broken powerlaw model 

The fit by means of the Broken powerlaw model does not provide a good description of the 
data, but still the variation of the indexes Fi and Fi between the states A and B can give some 
some clue on the physics. The indexes Fi and F2 are clearly different between spectra A and B, 
leading to a variation of the spectral break AF = F2 — Fi from 0.8 (spectrum B) to 1.1 (spectrum 
A), and the break energy Et, seems to be larger in state A than in state B, even if it cannot be well 
constrained because it lies in the gap between the data (between SPI and LAT). Spectrum C instead 
looks quite similar to spectrum B, except for the normalization which is lower. 

The indexes variation between the spectra suggests a modification of the electron distribution 
in the jet during the 7-ray flare. 
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Figure 2: 3C 273 multiwavelength spectra fitted with Broken powerlaw model (left panel) and Cutoff 
powerlaw+Log-parabola model (right panel). Solid lines are the best fit models. In the right panel the 
individual components Cutoff powerlaw and Log-parabola are plotted with dashed lines. Blue: spectrum A; 
Yellow: spectrum B; Green: spectrum C. See text for details. 



3.2 Log-parabola and Cutoff powerlaw+Log-parabola models 

The Log -parabola is the common model used to fit blazar spectra in /-ray energies [4]. The 
variation of the spectral slope of the Log-parabola model can represent much better the Fermi 
observations obtained for the spectrum A, when compared to a powerlaw with a constant slope. 

The bump in the high energy part of the spectrum is commonly considered to be due to Inverse 
Compton (IC) processes, could be Synchrotron Self Compton (SSC) or External Compton (EC), 
on the high energy electrons in the jet and the Log-parabola shape is able to model it. The shape 
may be given by a particle acceleration due to a stochastic component in the acceleration process 
[6]. 

The Log-parabola model does not provide a good description of spectrum A (Xred ^ ^)' ^^" 
cause the intrinsic curvature of Fermi data in spectrum A does not match with the X-ray data in a 
Log-parabola shape. Assuming that the Log-parabola is the model for the IC bump this suggests 
that there is an extra component in the spectrum. 

The presence of two components can be also inferred by looking at the lightcurves in Fig. 1 . 
The variability pattern of the two lightcurves at different energy bands is very different, suggesting 
that X and y-rays are not correlated and therefore they should be two distinct components or at 
least a single component with some variable parameter. In Soldi et al. 2008 [5] 3C 273 variability 
is discussed: they suggest two distinct physical components or a single component with varying 
spectral shape in the X-ray emission below and above ~ 20keV based on the lightcurve variability. 

We model the extra component in the spectrum as a Cutoff powerlaw at lower energies, for 
which we fixed the parameters in the fit. The adding of the Cutoff powerlaw allows us to im- 
prove the quality of the fit. The chosen fixed values of the Cutoff powerlaw parameters are indeed 
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Broken powerlaw 




Spectrum 


Ti 


Eh{MeV) 


Ti 


Nbp 


^red 


d.o.f. 




A 


1.57±0.13 


s.stli 


2.61 ±0.04 


(2.82±0.1)10-2 


1.442 


51 




B 


1.67±0.01 


l-6;ii 


2.41 ±0.85 


(4. 17 ±0.04) 10-2 


2.176 


43 




C 


1.62±0.01 


2.21^:^ 


2 43+°-^^ 


(0.91 ±0.03)10-2 


0.718 


16 






Log-parabola 






Spectrum 


a 


P 


N,p 


f^red 


d.o.f. 




A 


1.22±0.01 


0.108 ±0.002 


(1.54 ±0.04) 10-2 


2.233 


52 




B 


1.41 ±0.02 


(9.74;°-|)10-2 


(2.91 ±0.19)10-2 


0.965 


41 




C 


1.38±0.03 


(9.8i;°-^^)10-2 


(0.70 ±0.03) 10-2 


1.290 


17 






Cutoff powerlaw+Log-parabola 


Spectrum 


r 


EeikeV) 


Ncp 


a 


P 


Nip 


'^red 


d.o.f. 


A 


1.7 


150 


0.028 


76+°°'' 

^•'0-0.08 


0.16±0.01 


(2.0±0.4)10-3 


1.151 


49 


B 


1.7 


150 


0.028 


1 in+°°5 


n+ooi 

"•^-'-0.007 


(6.3 ±1.0)10-3 


1.132 


41 


C 


1.7 


150 


0.009 


nf, 1+0.20 


0.19 ±0.03 


(0.34;°i^)10-3 


0.098 


16 



Table 1: Table of fit parameters. The first column identifies the spectrum as explained in the text. Parameters, 
reduced X^ ^nd degree of freedom are reported for each spectrum and each model. Ni,p, Nip, Ncp are the 
normalizations of the respective models, corresponding to the photon flux at 1 keV. Errors are at the 90% 
confidence level (2.71 a) and the parameters without error are used as fixed parameters in the fit. 



quite arbitrary, because if we try to fit the whole set of 6 parameters of the Cutoff powerlaw+Log- 
parabola model there is some degeneration and it is not possible to find well constrained values, 
especially for the cutoff energy Ec which can span up to ~ 500 keV and still give a good fit. 

The F-test between the two models, Cutoff powerlaw and Cutoff powerlaw+Log-parabola 
on dataset A suggests that the second model is significantly better: the probability that Cutoff 
powerlaw+Log-parabola model is better by chance only is P/ = 3 • 10-^. 

The spectra B and C instead can be well fitted with a Log-parabola, because for these datasets 
the statistics at £ > 1 GeV (Fermi data in state B and EGRET data in state C) is poor, and the whole 
dataset is needed to estimate the curvature parameter j3. In this case is not possible to state if there 
is an extra component or not from the spectrum as in spectrum A. Using the Cutoff powerlaw+Log- 
parabola model with the same fixed parameters of spectrum A the fit in spectrum B and C is still 
good. 

The peak energy of the Log-parabola SED is given by Epeak = 10^2-«)/2p j^gy pj -pj^g peak 
energy of the three spectra seems to show a trend similar to that of Ey in the Broken powerlaw 
model, with a larger Epeak in spectrum A, and comparable values in spectra B and C. In any case 
the error on these parameters is quite broad so this trend may not be significant. 

We also tried to fit the three datasets A, B and C togheter with the Cutoff powerlaw+Log- 
parabola and the Log-parabola models, forcing the spectral parameters to be the same but leav- 
ing the normalization free for intercalibration. The different datasets can not be modeled by the 
same shape, proving that the difference between the datasets can not be due to the normalization 
only. The global fit of the two datasets A and B togheter has a x^ed — 13 using both the Cutoff 
powerlaw+Log-parabola and the Log-parabola model. 
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4. Conclusion 

A careful analysis is still in progress and more models can be tested. With respect to previ- 
ous work like [2] or [5] we have data in the y-ray domain up to ~ lOGeV from the Fermi-LAT 
spacecraft, allowing us to study the multiwavelength emission at higher energies with more detail. 

A possible and interesting scenario to interpret the data is the two components scenario, in 
which the y-ray emission modelled with the Log-parabola spectral law could be interpreted as a 
blazar-like emission due to IC processes (SSC or EC) in the jet, as suggested by the correlations 
between the hard X-ray and radio or optical emission claimed in [5]. 

The X-ray component we modelled with the Cutoff powerlaw could be a Seyfert-like compo- 
nent, due to thermal photon emitted by the accretion disk and scattered in the X-ray domain by IC 
scattering on a different electron population, possibly an hot corona close to the central black hole. 

More detailed analysis on the lightcurves and on optical and infrared data are needed to con- 
firm this scenario. 
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